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Abstract 

Objectives: Cervicovaginal fluid (CVF) can be considered as a potential source of biomarkers for diseases of the lower 
female reproductive tract. The fluid can easily be collected, thereby offering new opportunities such as the development of 
self tests. Our objective was to identify a CVF protein biomarker for cervical cancer or its precancerous state. 

Methods: A differential proteomics study was set up using CVF samples from healthy and precancerous women. Label-free 
spectral counting was applied to quantify protein abundances. 

Results: The proteome analysis revealed 16 candidate biomarkers of which alpha-actinin-4 {p = 0.001) and pyruvate kinase 
isozyme Ml /M2 (p = 0.014) were most promising. Verification of alphia-actinin-4 by ELISA (n = 28) showed that this candidate 
biomarker discriminated between samples from healthy and both low-risk and high-risk HPV-infected women (p = 0.009). 
Additional analysis of longitudinal samples {n = 29) showed that aiplia-actinin-4 levels correlated with virus persistence and 
clearing, with a discrimination of approximately 18 pg/ml. 

Conclusions: Our results show that CVF is an excellent source of protein biomarkers for detection of lower female genital 
tract pathologies and that alpha-actinin-4 derived from CVF is a promising candidate biomarker for the precancerous state 
of cervical cancer. Further studies regarding sensitivity and specificity of this biomarker will demonstrate its utility for 
improving current screening programs and/or its use for a cervical cancer self-diagnosis test. 
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Introduction 

The human papillomavirus (HPV) is responsible for virtually all 
cervical cancers [1]. Although more than 150 variants of this virus 
exist, only certain genotypes, such as HPV 16, 18, 33, 45 and 58 
are known as high-risk types (HR-HPV) [2] . Low-risk HPV types 
(LR-HPV), mainly HPV 6 and 1 1 , seldom cause genital tumors; 
however, they do cause condylomata acuminata (anogenital warts) 
[3]. Persistent HPV oncoprotein expression (E6/E7) in HPV 
infected epithelial basal cells deregulates cell division [4]. 
Overexpression of these viral genes causes the deregulation of 
cell proliferation, metabolism, apoptosis, differentiation and 
genomic instability, all of which may lead to consecutive stages 
of cervical intra-epithelial neoplasia (CINl, 2 and 3) or squamous 
intra-epithelial lesions (low-grade SIL and high-grade SIL). 
Approximately 85% of all sexually active persons will be exposed 



to HPV, and the majority of all HR-HPV infections are virion 
producing (which is limited in time) suggesting that in addition to 
the viral genotype, several cofactors are closely associated with the 
persistence of the viral oncoproteins and the transformation of the 
cervical mucosa to malignant tissue [5,6] . 

There are several reasons why the development of a diagnostic 
assay, based on cervicovaginal fluid (CVF), for early detection of 
cervical cancer would be beneficial. At first, current screening 
assays are not optimal. Because of the high prevalence of cervical 
cancer, females between the ages of 25 and 65 are frequently 
screened using Pap(anicolaou) smears, a test which is based on 
detection of morphological changes of cervical mucosal cells. 
Unfortunately, the low sensitivity of this test can result in cervical 
cancer diagnosed at a late stage [7]. Data on randomized trials 
show that HPV DNA screening is more sensitive than cytology, 
but specificity is lower, which inevitably will result in an 
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overtreatment because about 80% of the HPV positive patients 
spontaneously clear the virus [3]. Additionally, colposcopic 
examination is labor-intensive, difficult to automate and 

vulnerable to inter- and intra-observer variation [8]. The 
introduction of a new biomarker that helps improving sensitivity 
and specificity of current screening programs is therefore more 
than welcome. The second reason lies in the applicability of 
CVF for diagnosis. Because the fluid can easily be collected by 
the individual with the aid of special devices [9], a fast and 
simple test on the basis of a dipstick assay, performed by the 
woman herself, could be developed. Such a self test could be 
introduced e.g. in resource-limited populations where unequal 
burden of cervical cancer often occurs. Indeed, it is estimated 
that only 5% of women in low-resource countries are screened 
appropriately for cervical cancer [10]. Lack of healthcare 
infrastructure and financial cost are the main reasons why 
cytology-based programs are not implemented. Since alternative 
methods are currently investigated in these regions, use of a 
CVF-based self-diagnosis test could be considered. In addition, 
self tests could also be used in follow-up studies of vaccination 
programs since the effect of vac[:ination in sexually active 
women is uncertain [11]. Moreover, current vaccines only 
address neoplasia-inducing genotypes HPV 16 and —18 suggest- 
ing that under ideal circumstances, only a maximum of 70% of 
all cervical cancers can be prevented [7,12]. Thus at least for the 
coming decades, careful evaluation of \'a( ('inati()n programs is 
highly recommended and a simple self-diagnosis test could help 
a lot in meeting this question. 

In the past, several proteomics studies were performed to 
identify biomarkers for the (early) detection of cervical neoplasia. 
In different studies, precancerous and healthy cervical tissue was 
compared and Lee el al. studied the in viiro contribution of 
oncoprotein E7 to the development of cervical neoplasia using a 
HPV-infected cell Une [13-20]. Several markers that were Unked 
to dysregulation of the normal cell cycle, such as pi 6, Ki67 and 
cyclin E were characterized. These findings resulted in improved 
understanding of HPV-induced carcinogenesis but so far they all 
slio\\Td low sensitivity and/or specificity and therefore could not 
be used in the clinic. 

For biomarker discovery and validation, body fluids are better 
than tissue samples because the following advantages are inherent 
to most body fluids: (1) easy accessibility; (2) avoidance of sampling 
risks; and (3) multiple sampling potential [21]. Today, plasma is 
the most analyzed biological fluid in the search for biomarkers and 
it has already been used in cervical cancer research [22,23]. 
However, plasma comes into contact with nearly all organs of the 
body, making plasma biomarkers less specific because it is difficult 
to determine from which organ they originate [21]. 

CVF can be obtained by cervicovaginal lavages, vaginal swabs, 
tampons or with a device for self-sampling [9,24] and has the 
following benefits over plasma: (1) it possesses a lower fluid 
volume /organ-ratio compared to plasma, resulting in iiigher 
sensiti\'ity; and (2) it comes into contact with fewer organ systems, 
which most likely increases the specificity of the markers [21,25]. 
In this study, we focus on the identification of candidate 
biomarkers for cervical neoplasia by comparing the protein 
composition of individual CVF samples from healthy and 
precancerous (LSIL and HSIL) women. Our findings indicate 
that alpha-actinin-4 (ACTN4) is a promising candidate biomarker 
for the precancerous state of cervical cancer and that this CVF 
protein is very well suited for the development of a self-diagnosis 
test. 



Materials and Methods 

Study design and sample collection 

All patients agreed to participate by written consent and the 
study was approved by the ethical committee of the Antwerp 
university hospital (registration number: B30020 108372). Study- 
specific patient identification codes were assigned and transmitted 
in such a manner that patient confidentiality was preserved. 

AU samples in this prospective, blinded, cohort study were taken 
by the same gynaecologist (WAAT). In case of abnormal Pap 
smear results (due to HPV infection or as a false positive outcome), 
both healthy and precancerous patients are routinely subjected to 
a colposcopic examination, a procedure that includes rinsing the 
vagina with 5% acetic acid. After colposcopy this washing fluid 
(containing the cervicovaginal fluid) is normally discarded but was 
collected for proteomic analysis. Additionally, cervical cytology 
samples were collected during this colposcopic examination to 
determine the cytology and HPV status as previously described by 
using type-specific PCRs [26]. AU BD-SurePath liquid based 
cytology samples were sent to the pathology laboratory (RIATOL, 
Department of Molecular Diagnostics, Sonic Healthcare Benelux, 
Antwerp, Belgium) for processing. Based on the combination of 
colposcopic results, cytolog)- outcome and HPV genotyping, 
healthy and precancerous women could be identified on a reliable 
manner. We selected samples originating from 6 healthy (normal 
colposcopy/cytology and HPV negative; group A) women and 6 
precancerous (abnormal colposcopy/cytology and HR-HPV 
positive; group B) individuals (Tabk- 1). All samples were derived 
from postmenopausal (> three years after the last menses) women 
from similar age (59+/— 13 years). Patients who used medication 
were excluded from the study. Besides, patients who showed 
visible signs of gynaecological infections (e.g. gonorrhoea) or 
suffered a HIV-infection were excluded as well. Based on the 
cervicovaginal fluid samples, patients were screened for trichomo- 
nas vaginalis. During proteomic analyses, healthy and precancer- 
ous CVF samples were always analyzed in pairs to minimize 
technical variations of the LC-MS/MS platform. 

Proteomic analysis 

After collection, the cervicovaginal lavages (25-40 ml) were 
immediately transported on ice to the laboratory and stored at 
— 80°C. After centrifugation, the supernatant was concentrated by 
lyophUization to a final volume of approximately 200 (xL. Protein 
samples (1 mg/sample) were separated and fractionated in the first 
dimension on a reverse phase (RP) protein C4 HPLC column. 
Enzymatic digestion with tr)psin was performed and the resulting 
peptides of each fraction were separated in a second dimension on 
a RP-C18 micro-capillary HPLC system. Mass spectrometric 
analysis was performed using a MALDI-ToF/ToF instrument. 
Resulting MS/MS Spectra from each sample were screened 
against the human Swiss-Prot database (version: 57.1) using the 
MASCOT search engine. Analysis of the obtained datasets were 
performed as previously described [25]. More detailed information 
can be found in the Document SI and Figure SI. 

Enzyme linked immunosorbent assay 

To confirm the LC-MS/MS results, c'nz5'me linked immune 
sorbent assays (ELISA) were performed for alpha-aclinin-4 
(ACTN4) and pyruvate kinase isozyme M1/M2 (PKM2), according 
to the manufacturer's instructions. A human ACTN4 ELISA kit 
(Cusabio Biotech Co. LTD.) and a human PKM2 ELISA kit 
(USCNK Life Science Inc.) were used to quantify ACTN4 and 
PKM2 in cervical vaginal washes from healthy (n = 16) and HPV- 
infected (n=12) females (both pre- and post-menopausal). 
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Table 1. Patient information of samples used in the comprehensive proteomic experiments. 



Ssmpie 




G6notype(s) 


Virsl loBcl 


Colposcopy 


Cytology 


Unique Pept 


A.l 


57 


no HPV 




Normal 


Normal 


/ 


A.2 


69 


no HPV 




Normal 


Normal 


/ 


A.3 


52 


no HPV 




Normal 


Normal 


/ 


A.4 


46 


no HPV 




Normal 


Normal 


/ 


A C 
n.j 


00 


no nrv 




Normal 


Normal 


/ 


A.6 


62 


no HPV 




Normal 


Normal 


/ 


B.1 


61 


31/52 


335/2 


LSIL 


CIN1 


4 


B.2 


53 


16 


417 


HSIL 


CIN3 


2 


B.3 


47 


39 


2011 


LSIL 


CIN2 


3 


B.4 


72 


;s/56 


0.1322/1119 


LSIL 


CINl 


2 


B.5 


57 


31 


1148 


HSIL 


CIN3 


12 


B.6 


60 


31/52 


258/25 


LSIL 


CIN2 


6 



Samples Al-6 are originating from healthy Individuals {group A), while samples Bl-6 originates from precancerous patients {group B). HPV viral load is expressed as 
number of copies per cell. Colposcopy: Normal = no aberrant lesions detected, LSIL = low-grade squamous intra-epithelial lesions, HSIL - high-grade squamous intra- 
epithelial lesions, ASCUS = atypical squamous cells of undetermined significance. Cytology is listed according to the Cervical Intraepithelial Neoplasia {CIN) 
nomenclature. The numbers of unique peptides of ACTN4 are listed in the last column. 
doi:l 0.1 371 /journal.pone.Ol 06488.t001 



Additionally, several longitudinal samples (n = 29) were also tested 
for these two candidate biomarkers. 

Statistics 

To determine whether the difTerence in protein abundances was 
significant, normalized spectra abundance factors (NSAF values) 
were subjected to statistical testing. Because the selected proteins 
were identified in several samples, and therefore also contained 
several NSAF-values, a non-parametric Mann- Whitney U-test was 
used. Moreover, a chi-square test was performed, whereby the 
frequency of all identified proteins was analyzed to determine 
whether certain proteins were more or exclusively detected more 
under certain conditions. Unpaired Student's T-tests were 
performed to analyze the ELISA results. AU statistical tests were 
performed using the Statistical Package of Social Science (SPSS 
version 18). 

Results 

Differential study for the characterization of cervical 
cancer biomarkers 

To identify differentially abundant proteins, 12 samples were 
divided in two groups and were individually analyzed using an 
analytical proteomics platform (Table 1). Group A comprised of 6 
samples from healthy individuals, whereas group B consisted of 6 
samples from precancerous individuals. Proteomic analysis result- 
ed in 846 and 825 identifications for group A and group B, 
respectively, representing 371 (group A) and 341 (group B) non- 
redundant proteins (Figure 1). The mean overlap between all 
samples from group A was 54% and from group B was 53% 
(Table SI). 

A semi-quantitative analysis was also performed based on the 
label-free NSAF method (Table S2). For both groups, the NSAF 
values of the proteins present in at least 5 of the 6 samples were 
mutually compared. The variation in protein abundance between 
samples was expressed as the coefficient of variation (CV). A mean 
CV value of 60% (lowest: 19%; highest: 134%) was obtained for 
group A (healthy), while a mean CV of 52% (lowest: 14%; highest: 
99%) was noted for group B. 



Identification of differentially abundant proteins 

To identify candidate biomarkers that correlate with the 
precancerous state of cervical cancer, healthy and precancerous 
CVF proteomes were compared. A qualitative analysis was 
performed between the 371 identified proteins for the healthy 
group and the 341 identifications for the precancerous group, 
which resulted in 238 or 67% overlapping proteins. To identify 
proteins with a statistically significant difference of occurrence, all 
identifications were subjected to statistical testing using a Pearson 
chi-square test. Proteins with a p value ^0.05 are listed (Table 2). 
From these calculations, ACTN4 was identified in all samples 
from precancerous patients but not in the samples from healthy 
patients. Moreover, ACTN4 was primarily identified based on 
several unique peptides with a mean MASCOT protein score of 
162 (Table S3). The difference in the abundance of this protein 
was further verified by ELISA (see below). 

In addition, a semi-quantitative comparison was performed to 
determine whether certain proteins that were identified in both the 
healthy and the precancerous groups were differentially expressed. 
The NSAF-values of corresponding proteins were statistically 
tested using a Mann-Whitaey U-test. Only proteins with a 




Figure 1. Venn diagram of the healthy and precancerous 
proteomes. For the healthy proteome 371 unique proteins were 
identified while 341 unique proteomes for the precancerous proteome. 
An overlapping set of 238 proteins were found in both proteomes. 
doi:1 0.1 371/journal.pone.01 06488.g001 
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Table 2. Proteins with a statistical significant difference in frequency between the 


healthy and preca 


ncerous group. 




Acc.No 


Protein 


p-value 


# Healthy 


# Precanc 


043707 


Alpha-actinin-4 


0.001 


0 


6 


P00558 


Phosphoglycerate kinase 1 


0.014 


2 


6 


P02774 


Vitamin D-binding protein precursor 


0.046 


6 


3 


P08603 


Complement factor H 


0.046 


3 


0 


P29373 


Cellular retinoic acid-binding protein 2 


0,046 


0 


3 


P29508 


Squamous cell carcinoma antigen 1 {SCCA-1); Serpin B3 


0.046 


3 


6 


P43490 


Nicotinamide phosphoribosyltransferase 


0.046 


0 


3 


P61158 


Actin-related protein 3 


0,046 


0 


3 


P62258 


14-3-3 protein epsilon 


0,014 


0 


4 


P62988 


Ubiquitin 


0.014 


4 


0 


Q5VTE0 


Putative elongation factor 1 -alpha-like 3 


0,046 


3 


6 


Q9UIV8 


Serpin B13 


0.046 


0 


3 


Acc No represents the accession number of the identified proteins. Statistical significance 
show the number of samples that contain the corresponding protein. 
doi:l 0.1 371 /journal.pone.Ol 06488.t002 


is listed in column named p-value. Columns "j^ Healthy" and "# Precanc" 



Statistical significance of p^O.05 are listed (Table 3). This 
calculation resulted in four proteins that were characterized by a 
significant difference in abundance between both conditions. Only 
one protein, haptoglobin, showed a downregulation in the 
precancerous condition, whereas three other proteins, ATP 
synthase subunit beta, annexin A2 and PKM2, were upregulated. 
PKM2 and annexin A2 showed the lowest p value (0.01 and 0.03, 
respectively) and were detected in 5 of 6 samples from the 
precancerous group. Because highly variable proteins are not 
suitable candidate biomarkers, the variation in the abundances of 
these proteins was calculated. With CV values of 40% and 29% 
for the healthy and precancerous group for PKM2, respectively, 
this protein is more suitable as a candidate biomarker compared to 
annexin A2, which had values of 77% and 48%, respectively. 

Verification of the candidate biomarkers 

An ELISA was performed to verify the differential profiles of 
ACTN4 and PKM2. Healthy (n= 16) and HPV-infected samples 
(both HR-HPV oncogenic (n = 8) and LR-HPV non-oncogenic 
samples (n = 4)), which were not included in the previous 
experiment, were tested using commercially available kits 
(Table 4). 

Levels of ACTN4 were significantly higher in samples from 
HR-HPV infected individuals compared to samples from healthy, 
non-HR-HPV infected patients (p = 0.023), which confirmed our 
LC-MS/MS results. Remarkably, sample one, originating from a 



non-HR-HPV infected patient, appears to be an outlier because of 
the high ACTN4 concentration of 85.6 pg/mL (see discussion). 
Additionally, the LR-HPV infected patients had higher ACTN4 
concentrations compared to the non-HPV-infected patients 
(p = 0.052), but somewhat lower values compared to the HR- 
HPV infected patients (p = 0.114). Comparison of the ACTN4 
levels in HPV-infected (both HR- and LR-infections) versus non- 
infected samples (healthy patients) showed a significant difference, 
with a p value of 0.009 (Figure 2). The highest ACTN4 
concentration for the healthy group (except for the outlier) was 
17.3 pg/ ml, whereas the lowest concentration for the HR- infected 
group was 30.6 pg/ ml. The lowest concentration of ACTN4 in 
the LR- infected group was 21.3 pg/ml, which is higher than the 
highest value of the healthy group. Based on these observations, a 
threshold of 1 8 pg/ml can be set as the highest value of ACTN4 
for the non-infected condition. No linear correlation was found 
between the viral load and the absolute concentration of ACTN4, 
which presumes that ACTN4 is not just a marker for HPV 
infection but rather for precancerous stage. 

In addition, the most promising semi-quantitative candidate, 
PKM2, was verified with a commercial ELISA kit. For this 
analysis, the same 28 samples were used; however, the LC-MS/ 
MS results were not confirmed (data not shown). 

We then analyzed the ACTN4 concentrations for an additional 
29 longitudinal samples (i.e., samples from the same patient at 
several time points) from 9 patients. For each patient, a minimum 



Table 3. Proteins with a statistical significant difference in abundance between the healthy and precancerous group. 



Acc.No Protein p-value Mean Healthy Mean Precanc. Ratio 



P00738 


Haptoglobin 


0.03 


2,630 (3) 


1,227 (5) 


47% 


P06576 


ATP synthase subunit beta, mitochondrial 


0.05 


0,816 (3) 


4,825 (3) 


591% 


P07355 


Annexin A2 


0.03 


3,831 (6) 


11,266 (5) 


294% 


P14618 


Pyruvate kinase isozymes M1/M2 


0.01 


0,4S2 (4) 


1,401 (5) 


310% 



Columns "Mean Healthy" and "Mean Precanc." show the mean normalized spectral abundance factors (NSAF-value) while (x) indicates the number of samples that 
contain this ID. Ratio is expressed as precancerous/healthy. "Acc No" represents the accession number of the identified proteins and "p-value" shows the statistical 
significance. 

doi:1 0.1 371/journal.pone.01 06488.t003 
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Table 4. Patient information of samples used for the validation of ACTN4 by ELISA. 



Condition 


Age 


Genotype(s) 


Viral load 


Colposcopy 


Cytology 


Cone ACTN4 


Healthy 1 


54 


no HPV 




Normal 


Normal 


85.6 


Healthy 2 


53 


no HPV 




Normal 


Normal 


1.1 


Healthy 3 


41 


no HPV 




Normal 


Normal 


12.2 


Healthy 4 


21 


no HPV 




Normal 


Normal 


3.7 


Healthy 5 


50 


no HPV 




Normal 


Normal 


14.4 


Healthy 6 


52 


no HPV 




Normal 


Normal 


1.0 


Healthy 7 


59 


no HPV 




Normal 


Normal 


6.2 


Healthy 8 


34 


no HPV 




Normal 


Normal 


7.2 


Healthy 9 


36 


no HPV 




Normal 


Normal 


-0.2 


Healthy 10 


31 


no HPV 




Normal 


Normal 


10.6 


Healthy 1 1 


51 


no HPV 




Normal 


Normal 


0.2 


Healthy 12 


30 


no HPV 




Normal 


Normal 


17.3 


Healthy 13 


28 


no HPV 




Normal 


Normal 


15.0 


Healthy 14 


37 


no HPV 




Normal 


Normal 


6.9 


Healthy 15 


38 


no HPV 




Normal 


Normal 


-13.0 


Healthy 16 


62 


no HPV 




Normal 


Normal 


16.3 


Low Risk 1 


33 


6 


96249 


ASCUS 


Normal 


25.4 


Low Risk 2 


37 


11 


51740 


ASCUS 


CINl 


21.3 


Low Risk 3 


42 


6 


1 


Normal 


Normal 


41.3 


Low Risk 4 


49 


6 


0.05 


Normal 


Normal 


59.8 


High Risk 1 


35 


16/39 


7729/1661 


ASCUS 


CIN3 


34.1 


High Risk 2 


51 


52 


129 


LSIL 


CINl 


227.8 


High Risk 3 


20 


16/31/52/66 


0.02/11/12/31 


LSIL 


CINl 


65.2 


High Risk 4 


32 


16/31/39/52/66 


288/0.2/4/7416/179 


LSIL 


CINl 


44.4 


High Risk 5 


31 


16/58 


11571/253 


HSIL 


C1N2 


89.9 


High Risk 6 


32 


16/53/58/59 


9126/79/2733/1510 


LSIL 


CINl 


140.8 


High Risk 7 


28 


35 


5159 


LSIL 


CIN2 


30.6 


High Risk 8 


33 


31 


1696 


HSIL 


CINl 


30.9 



HPV viral load was expressed as number of copies per cell. Colposcopy: Normal = no aberrant lesions detected, LSIL = low-grade squamous intra-epithelial lesions, 
HSIL = high-grade squamous intra-epithelial lesions, ASCUS = atypical squamous cells of undertimined significance. Cytology results were listed according to the 
Cervical Intraepithelial Neoplasia (CIN) nomenclature. The concentration of ACTN4 was described in pg/ml. 
doi:l 0.1 371 /journal.pone.Ol 06488.t004 



of three longitudinal samples was evaluated (Figure 3). From tiiese 
nine patients, five patients cleared the virus, two patients had 
persistent HPV infection, one patient acquired a new HPV 
infection and one patient was healthy. Although the numbers of 
patients per group are small, these samples can help us to confirm 
our findings. An ACTN4 ELISA was performed on several of the 
samples to identify correlations between the virus profiles and 
ACTN4 levels within one patient. However, in addition to the 
virus titer, we also expected the genotype and time of infection to 
influence the ACTN4 expression because different virus types 
have different carcinogenic capacities and cervical cancer usually 
develops after a persistent HR-HPV infection [2] . Recendy, it was 
shown that the development of cervical precancer (cervical 
intraepithelial neoplasia of grade three; CIN3) is preceded by a 
steady increase in the viral load of a given HR-HPV type, whereas 
a rapid exponentially increasing load is generally cleared within 6 
to 18 months and is usually associated with low-grade cytological 
abnormalities [27]. With these factors in consideration, the 
following two trends were observed from the experiments: (1) 
patients who had an early-stage infection (patient 1) or cleared the 
virus (patients 2, 6, 7 and 9) had increasing or descending ACTN4 



levels, respectively, and patients with no or low infections (patient 
.5) or continuous infections (patients 4 and 10) had stable low or 
high ACTN4 levels, respectively; and (2) for a persistent infection 
(patients 1, 4 and 10), the ACTN4 accumulated. 

Discussion 

Characterizing the CVF proteome of healthy and 
precancerous patients 

In this study the CVF proteome from 12 (6 healthy and 6 
precancerous) individual samples was characterized to identify 
candidate protein biomarkers for the detection of the precancerous 
state of cervical cancer. The analysis of these 12 samples was 
performed at both the qualitative (present/absent) and semi- 
quantitative level. The number of protein identifications for each 
of the samples is very consistent with previously performed 
proteomic studies on CVF [28]. 

The determination of the amount of overlapping proteins 
between the samples within one group showed significant 
diflFerences. Only 10% of all proteins were observed in each 
individual sample from both groups. This group consists primarily 
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Figure 2. ACTN4 levels between different study groups. Concentration of ACTN4 for the healthy individuals (n = 16), the HR-HPV infected 
patients (n = 8) and the LR-HPV infected patients (n = 4). Statistical testing of ACNT-4 levels between healthy and HR-HPV infected patients resulted in 
a p-value of 0.023 and between healthy versus LR-HPV in a p-value of 0.052. When the HR- and LR-HPV infected patients were grouped as HPV 
infected patients, a p-value of 0.009 was found between non-infected (healthy patients) and infected groups. 
doi:1 0.1 371 /journal.pone.01 06488.g002 



of proteins that are characteristic of tlie CVF proteome, based on 
their localization and biological function. We therefore call these 
proteins "characteristic proteins". However, a significant number 
of identifications, such as intracellular proteins, are less charac- 
teristic of the CVF proteome. This finding can be explained by 
processes such as the disruption of the epithelial cell layer during 
sampling and the shedding of dead epithelial cells. We therefore 
group these identifications as "non-characteristic" proteins. Based 
on these findings, which are common to almost every compre- 
hensive CVF proteomics study, one can make a distinction 
between a "core proteome" of frequently identified characteristic 
CVF proteins and a "variable proteome" that contains the non- 
characteristic CVF proteins [28] . 

Compared to the mean CV values of the common identifica- 
tions between three technical replicates (25%) [25], the mean CV 
values of protein abundance in the healthy and the precancerous 
samples (60% and 52%, respectively) were significantly higher. 
This result was expected because, in contrast to the technical 
replicates, the samples were not identical. In all samples, proteins 
with high inter-individual differences were observed (up to 134%), 
whereas the measured mean CV value of the "characteristic 
proteins" was 38%. These results are similar to other proteomic 
studies of body fluids. For example, Yamakawa el al. [29] found a 
mean CV value of 63% for the seminal plasma proteome, and 



according to Anderson et al. [30] , the mean CV value of a series of 
plasma proteins was approximately 45%. 

Identification of differentially abundant proteins 

To examine differential abundant proteins, a label-free spectral 
counting method was performed. As this method allows only for a 
semi-quantitative analysis, outstanding effects were first investi- 
gated as these usually give the most rehable results. Therefore, we 
first pursued proteins for which their presence or absence in the 
precancerous group was characteristic for a given condition. 
Statistical testing (chi-square test) resulted in 12 proteins with a 
significant difference of incidence (p<0.05), of which 6 proteins 
were exclusively identified in precancerous samples (Table 2). The 
most significant protein was ACTN4 (Acc No: 043707), which 
was identified in all precancerous samples but none of the healthy 
samples (p = 0.001). Notably, SCCA-1, often used as plasma 
biomarker for cervical cancer [31], was also identified, albeit with 
a higher p value (p - 0.046). 

ACTN4 is a critical component of the cytoskeleton and is 
involved in the formation of cell-cell adhesions [32]. It has been 
shown that the inhibition of certain membrane bound adhesion 
molecules caused by a malignant transformation results in a 
reduction of cell-cell adhesion strength and upregulation of 
ACTN4 [33]. Based on a knockdown study of ACTN4 using 
breast cancer cells, Khurana el al. demonstrated that ACTN4 was 
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Figure 3. Evolution of ACTN4 levels and HPV infection on longitudinal samples. Comparison of ACTN4 levels with genotype and viral load 
in longitudinal samples from 9 different patients. Similar genotypes are illustrated in the same colours. At the left side the viral load (in number of 
copies/cell) is shown on a logarithmical scale while at the right side the concentration of ACTN4 is expressed in pg/ml. The y-axe indicates the time 
(expressed in days) while the patients were followed. A black dotted line illustrates the ACTN4 threshold of 17.3 pg/ml. ACTN4 levels were 
determined on the time points when a CVF sample was available. 
doi:1 0.1 371/journal.pone.01 06488.g003 



involved in the control mechanism of cell growth [34] . Moreover, 
the cytoplasmic accumulation of ACTN4 was associated with 
several malignancies and was linked with a bad prognosis [35], 
and an accumulative genomic gain of ACTN4 was associated with 
the formation of ovarian adenocarcinomas [36]. Also, RNAi- 
mediated downregulation of ACTN4 showed an association 
between the over expression of ACTN4 and an aggressive 
phenotype of oral squamous cell carcinoma [37]. Recently, the 
protein was found to be associated with several factors with 
varying functions, suggesting other roles for this protein, apart 
from actin regulation [38]. Indeed, Khurana et al. identified a 
LXXLL motif that was responsible for interaction of ACTN4 with 
the estrogen receptor and co-activators [39]. Recent studies in 
HPV transgenic mouse models provide evidence that estrogen and 
its nuclear receptor promote cervical cancer in combination with 
HPV oncogenes [40]. It is therefore tempting to speculate that 
ACTN4 plays a role in this important effect. 

Based on the semi-quantitative NSAF-method, proteins were 
selected for their differential expression between the two condi- 
tions. Statistical testing (Mann-Whitney U) resulted in four 
proteins with a significant (p<0.05) difference in abundance 
(Table 3). Although these four proteins show a statistical 
significant difference during this discovery phase, a verification 
step is needed to confirm these results. When both the p-value and 
the frequency of presence were taken into account, annexin A2 
and PKM2 were most promising (p<0.03) because both were 
identified in at least 9 of 12 individual samples, with a three-fold 
upregulation for the precancerous condition (Table 3). Despite the 
identification of annexin A2 in nearly every sample (11/12 



samples) with an acceptable difference in abundance (p = 0.03), the 
inter- individual variations (CV 77%) for the healthy group were 
much higher compared to PKM2 (35%). Therefore, PKM2 (Acc 
No: PI 46 18) was selected for verification by ELISA. 

PKM2 is a glycolytic enzyme that contributes to the energy 
supply of the cell. This protein has two isoforms, an embryological 
isoform (M2) with elevated enzymatic activity and an adult isoform 
(Ml). Recentiy, it was demonstrated that malignant transforma- 
tions of cells by HR-HPV oncogenes E6 and E7 induced a PKM 
switch from Ml to M2 [41]. This switch causes a shift from normal 
cellular metabolism to elevated glycolysis, which is beneficial for 
the growth and neoplasm of tumor cells [42]. Unfortunately, we 
were not able to confirm the difference in PKM2 concentrations 
between CVFs from healthy and precancerous individuals, 
although two ELISA kits from different suppliers were used. It is 
possible that the semi-quantitative spectral counting method may 
be prone to too much variation. Alternatively, we have previously 
observed that the CVF matrix is not optimal for some ELISA 
configurations and this may be the case for the two PKM2 ELISA 
assays. 

Verification of alpha-actinin-4 levels 

/Upha-actinin-4, the most promising candidate biomarker, was 
verified on 57 additional CVF samples (28 individual samples and 
29 longitudinal samples) by ELISA, a number that certainly fulfills 
the requirements for the discovery /verification phase of protein 
biomarkers [43]. The ELISA results showed that the concentra- 
tion of ACTN4 is elevated in the HPV infected patients compared 
to the non-infected patients, thereby confirming the LC-MS/MS 
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results. However, infection with LR-HPV (HPV-6 and HPV-U) 
resulted in higher ACTN4 concentrations in the CVF. Although 
the concentrations of ACTN4 measured in the CVFs from LR- 
HPV-infected individuals are higher than those from non-infected 
individuals, no significant difference in the abundance level was 
observed (p = 0.052). Because these findings are based on the 
analysis of samples from only four single LR-HPV-infected 
patients, more of this type of sample must be evaluated before 
conclusions can be drawn. One non-infected patient showed an 
outlier value of 85.6 pg/ml of ACTN4. This patient did not have a 
known history of HR-HPV infection, but the clinical file showed 
ischemic heart disease and a prevalence of breast cancer. It is 
therefore possible that this elevated concentration of ACTN4 is 
caused by the medical treatment or an LR-HPV infection other 
than HPV-6 or HPV- 1 1 . Therefore, extensive validation of this 
marker on a larger number of samples, including those from 
patients suffering from other diseases and infections (HIV, herpes, 
bacterial vaginosis. Chlamydia etc.), is the next step to determine 
whether ACTN4 is suitable for diagnosis purposes. In these 
studies, one could set a threshold at 18 pg/ml to discriminate 
between samples from healthy versus HPV-infected individuals. 

Although ACTN4 levels did not numerically correlate with the 
viral loads, an ascending and descending trend could clearly be 
observed in patients who either started an infection or cleared the 
virus. We also observed that infections over longer periods caused 
accumulation of the marker. A high level of ACTN4 may 
therefore be an indication of a persistent oncogenic HPV infection 
and an increased chance of developing cervical cancer [6,44]. It 
was generally observed in longitudinal samples that total viral 
loads of approximately 1 00 particles/ cell caused ACTN4 levels 
above the threshold value of 18 pg/ml. An exception on this rule 
was observed in patient 9, in which a viral load of 5,159 HPV35 
particles/cell resulted in only 12.7 pg/ml ACTN4 (Figure 3). 
Whether this false negative result can be explained by the very low 
oncogenic potential of HPV35 must be investigated [2]. 
Conversely, two patients infected with very low doses of HPV6 
had exceptionally high levels of ACTN4 (Table 4). These false 
positives could be the result of an exceptionally high ACTN4- 
inducing capacity by HPV6 or of an additional viral, fungal, or 
bacterial infection that also caused an increase of ACTN4 levels in 
the CVF. One must be aware however, that release of ACTN4 is 
possibly a consequence of an immune-induced lysis of a 
precancerous lesions, since HPV virions are shed into the 
environment in the absence of lysis or necrosis [45] . This suggests 
that HPV virus production and ACTN4 release do not necessarily 
occur simultaneously, although they usually do. In the two above 
mentioned patients however, virus production may have stopped 
while precancerous lesions or benign tumors continue to be target 
for the crippled immune system, restdting in a concomitant release 
of ACTN4. Further validation steps, on a larger amount of 
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